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Pricing American and Asian Options 

Pat Muldowney 


Abstract 

An analytic method for pricing American call options is provided; fol¬ 
lowed by an empirical method for pricing Asian call options. The method¬ 
ology is the pricing theory of M, which is also the source of the notation. 

Suppose a share has market value z{s) at time s, and suppose an option 
contract, between two parties, is entered into at time s = 0, so that one of the 
parties to the contract acquires the right, but not the obligation, to purchase 
the share, from the other party, for a fixed price k at any time t, 0 < t < r. 
The option expires at time r, so there is no right to purchase the share at times 
t > T. These are the features of an American call option. 

An American call option can be exercised at any time during the term or 
lifetime of the contract. This differs from a European call option which can only 
be exercised on the date of termination of the contract. 

Here are some questions about American call options. What is the initial 
value (or purchase price) r/;(0) of the option contract at time 0? At any time t, 
0 < f < T, what is the market value w{t) of the (unexercised) call option? 

If time t > T, then w(t) = 0 since the option can no longer be exercised. At 
time t = T, the American call (if not previously exercised) has value 

w{t) = max{z(r) — k, 0}, 

the same as the corresponding European call option. 

At time t < t, the American call option and the corresponding European 
call option have non-zero valu^ even if z{t) < k. The difference between the 
two is that, for t < t, the European call cannot be exercised even if z(t) > k. 

If z(0) is the observed value of the underlying share at the time of writing 
the contract, the purchase price w(0) and the exercise price k of this contract 
must satisfy 

w(0) + K > z(0). (1) 

The reason for this is that, if w(0) < z(0) — k, the option can be exercised 
immediately, the share can be purchased for the amount k, and then re-sold for 
amount 2:(0), realizing a riskless profit of z(0) — k — i(;(0). 

For t < T the American call option will not be exercised at time t if z{t) < k. 
What if z{t) > At? In that case, and disregarding the outlay required to purchase 
the option, the holder of the option will obtain immediate profit of amount 

^If K is excessively large the option may be effectively worthless. See Figure below. 
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z(t) — nhy purchasing the share for the option exercise price k and immediately 
re-selling the share for the market price z{t). If t = r that is what is expected to 
happen, since there are no further opportunities to purchase the share at option 
exercise price k. 

But if t < T, and even if z{t) > k, then, depending on various factors, the 
holder of the American call may decide not to exercise in case they expect the 
underlying share price to rise further in the time interval ]t,T]. If that were to 
happen, then the profit (or payojf) z(s) — k to be obtained by exercising the 
option at the later date s > t would be greater than what would be obtained 
by exercising at time t. 

This introduces a new element of unpredictability. If time t is the “present”, 
the share price z{t) has a known or observed value. Call this value Kt- Even if 
Kt > K, the option holder, in anticipation of further increase in the price of the 
underlying share, may decline to exercise the option at time t. 

On the face of it, in advance of exercising, the optimal time t at which to 
consider exercising the American call is the time t for which the share price z(t) 
is maximal, 0 < t < r; so x{t) > x(s) for s ^ t. 

However, the time value of money has to be factored into this argument. The 
time value of money is measured by the risk-free interest rate p. The following 
formulation takes account of time value of any potential option payoff. 

If r > t > s and 2 :(s) > z{t) > n, it is more advantageous to exercise the 
American call at the earlier time s if 

(z(s) - k) > z(i) - K. 

At time s = 0 the optimaj^ exercise time c is unpredictable. It is possible that 
the holder of an American call will exercise the option in advance of the optimal 
time or will do so after the optimal point has passed. Or will not exercise at 
all. 

Assuming “equal likelihood”, each of these possibilities is just as likely as 
the other. On this reasoning, the exercise time (or date) t of an American call 
option is a random variable whose expected value is c. 

To explore further, assume that t < r is the “present” so the share price 
z{t) = Kt is known. Suppose that, even if Kt > k, the option holder, in anticipa¬ 
tion of further increase in the price of the underlying share, declines to exercise 
the option at time t. Assume that the likelihood of the option holder making 
this decision is equal to the probability that z(s) > Kt for some s G ]t, r]. That 
is, 

2:(s) — K > z{t) — K = Kt — K. 

In reality it is necessary to factor in the time value of money here. The option 
holder will decline to exercise if they expect that the discounted value of the 
term on the left hand side of the following inequality exceeds the right hand 
side: 

(z(s) — k) > z{t) — K = Kt — K. 

^Likewise, at times s > 0 the optimal exercise time <; is unknown, s < <; < r. 
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If t is given, and if s is the “present”, with s < t, similar considerations apply 
to any decision to exercise the American call during the period ]s,t]. 

This argument implies certain assumptions about the motivation and be¬ 
haviour of option holders/buyers as a body, in regard to future prices of an 
underlying share. But the future market value of the underlying share is said to 
be “correctly” established by similar motivation and behaviour of those buying 
and selling the share. So there seems to be some justification for such assump¬ 
tions. 

Take time 0 to be the “present”. Suppose s is a given “future” time, 0 < 
s < r; so it is not known whether z(s) is going to be less than or greater than 
the exercise price k. Consider an option exercise decision at the future time s. 

For any exercise date s (0 < s < r) denote the discounted “difference” 
function for the option by 

p(s,z(s)) := (z(s) - (2) 

If time t > 0 is the present and 0 < t < s < r, then the relevant difference 
function is discounted from time s to time t: 

Pt{s, z(s)) := (z(s) - k) (3) 

Take time 0 to be the present, so option exercise decisions are based on (^. The 
option is exercisable (or “in the money”) if z(s) > k for some s g]0,t]; and it 
is optimally exercisable at time s = s{zt) if, for some other time s' before or 
after s, the payoff, discounted from time s' to time 0, is not greater than the 
payoff at time s(zt) discounted to time 0. The optimal exercise time s{zt), if 
exercis^is feasible, is variable for different histories zt- 

The graph or history zt is a potential occurrence of a joint-basic observable 
Zt if there is a distribution function F for whictj^ 

Zt zz Zt [R+ , F"] . 

In that case, for 0 < s < r the discounted payoff function values p(s, zt) are 
potential occurrences of a contingent observable 

pis, Zt) zi p(s, zt) [R+, F"] , 

whose expected value is 

<?(s) := E [p (s, Zt)] = / p{s, ZT)FiI[N]). (4) 

JRI 

If the exercise date s is selected, then s is fixed in this integration. Therefore 
the point-integrand n is a cylinder function, so the integral reduces to a one¬ 
dimensional integraP on R_j^ = 

^Exercise may never be feasible if the exercise price k is set too high. See Figure 4 below. 
^It is sometimes assumed, incorrectly, that F is the geometric Brownian distribution func¬ 
tion Q. 

^If F is the geometric Brownian distribution function then the integral is the one 
given in line 5 of m page 442, but with natural process growth rate /i in place of the risk-free 
growth rate p, and with p(s, zt) in place of maxfzir — /t, 0}. 
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Choose s so that q{s) is maximizecj^ and denote the chosen s by Then 
0 < <; < T, and, if time 0 is the present, then, at time 0, it is predicted that the 
American call option will be exercised at time 

The question also arises as to whether, if the “present” is some time t later 
than 0, the predicted exercise date will be different. 

The object of this article is to consider ways of finding a purchase price 
w{t),= Wt, for a call option at any time t, 0 < t < r; in particular the price 
■u;(0),= wq, payable by the initial purchaser of the option. Consider the latter 
problem first, so t = 0 is the present. 

For a European call option the exercise date is r, and the initial value of the 
option is obtained by means of the risk-neutral pricing formula (8.46) of [H] . 
page 440: 

e-"" / fr{zT)F{I[N]), 

Jul 

where F is the risk-neutral martingale distribution function discussed in Sections 
8.14 and 8.15, pages 436-440 of [T3], and = max{z(r) — /v,0}. 

For the corresponding American call option, the predicted exercise date is 
<j, which may be less than r. But apart from that, the pricing argument is the 
same and gives initial value of the American call as 


/ h{zT)F{I[N]), 

JTlI 


where f<;{zT) = max{z((r),0}. The integral on can be replaced by one on 
Either way, since the integrand is cylindrical the integral reduces to a 
one-dimensional integral on R+, = R^^^. 

When it is assumecQthat the underlying asset price process Zt is geometric 
Brownian with F = , then Section 8.16 ([14] pages 440-444) shows that 

F = , where a is the volatility, /i is the growth rate of the process, and p is 

the risk-free interest rate. In that case the initial value of the American call is 
given by (8.53) of [Tl] page 442, with r replaced by ■?. 

How can the optimal exercise time be established? If the process distri¬ 
bution function F is an expression such as then the expected discounted 
difference q{s) of has a particular form whose maximum value can be esti¬ 
mated; so the corresponding s, = can also be estimated. 

On the other hand, suppose F (or, at least, approximate values of F) is 
determined empirically, as indicated in m Section 9.11 page 488. Then it may 
be possible to calculate empirically the expected discounted difference q{s) of 
using different values of s, and use the results to estimate <^, the exercise 
date which gives the largest expected discounted difference. 

As an exercise in determining optimal exercise time q, take F = , which 

is the traditional (but generally incorrect) assumption. Then, using Lemma 33 

®Set aside, for the moment, the question of how to make this choice. 

^This assumption is generally incorrect. See m Section 9.9, pages 479-485. 
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(pages 309-310 of [13]) with change of variable z(s) = e“, Q becomes 


q(s)= [ e-^z{s)-n)g^^{I[N]y, (5) 

JrJ0,sJ 

SO 


<l{s) 




With /i = 0, cr = 0.5, T = 2, and with exercise price k = 0.1, 2, 1, and 5, 
Figures 1, 2, 3, and 4, respectively illustrate the expected values q{s) of the 
corresponding options if they are exercised at times s for 0 < s < 2. 

These graphs indicate that if exercise price k is very low (Figure 1), the 
situation can be compared to that of Q, and it is best to exercise the option 
immediately. In other words, at time t = 0 essentially riskless profit is available 
“here and now”; so why wait around for potential increases in the underlying 
share price? After all, the latter could actually decrease. 

If K is a bit larger (Figure 2), at time t = 0 it is safest to plan on exercising 
the option at or neaJj the contract termination date s = r. Because, even if 
(as in the cases portrayed here) the growth rate /i of the underlying share price 
is zero, there is a “structural” mean growth rate of due to volatility. So 
the attitude of the investor in American call options at time t = 0 will learj^ 
towards late exercise. 

For intermediate values of k (Figure 3) the optimal exercise time ^ « 1 is 
intermediate between 0 and 2. And for very large k (Figure 4), the long position 
in the option contract appears unattractive—potential buyers may not be too 
interested in the long position in an option with such a high exercise price, and 
will not pay very much for it; so it can only be sold for some very small amount 
of money. 

Figure 4 does not imply that an option with these parameters can never 
be in the money. The function q{s) yields expected rather than actual values. 
But the diagram suggests that the option should be immediately exercised if 
and when it achieves in-the-money status, without waiting around for further 
increases in underlying share price. Secondly, at time t = 0 the purchase price 
of such an option should be based on exercise at time s = t, not earlier, in order 
to allow for any possible growth in underlying share price. 

The optimal exercise date can be estimated by methods such as the one 
demonstrated in Figures 1 to 4. Once <: is established, the Black-Scholes pricing 
formulae (8.53) of [14] page 442 can be used, with c replacing r; giving the 

® Figure]^ gives a more detailed graph of g(s) for this case. It shows that the maximal value 
of g(s) occurs at about s = 1.83. 

®Some commentary on American call options suggests that there is no merit at all in early 
exercise. 
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Figure 1: /c = 0.1 Figure 4: k = 1 



Figure 2: k = 1.1 Figure 5: k = 5 



Figure 3: k = 1.1, i = 0, t = 2 Figure 6: k = 1.1, t = 1.5, r = 2 
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initial {t — 0) value of the American call option as 


Wo = 


<7^ 


— ne 


-P?N' 


0,1 




(6) 


Now suppose the option, created at time s = 0, is unexercised at time t, 0 < 
t < T. Suppose the holder of the option wishes to sell it at time t. What price 
should be charged for it? In other words, what is the value w(t) (or Wt)? 

Can it be expected that, at time t later than 0, the expected exercise date 
for the as yet unexercised option is the exercise date obtained by maximising 
q{s) in (U? If that were the case, then formula (8.54) of [T3] page 442 could be 
used for wt, with c replacing t; so Wt would be 


ztN' 


0.1 ( In t + (p + (<^ - t) \ 1 / \n^ + {p- - t) 


-1 


cry/<i - t 


( 7 ) 

The trouble with this argument is that is not necessarily the actual exercise 
date. It is a date estimated at time s = 0 for the purpose of pricing the option 
at that time. The option holder may choose to exercise at some time s earlier 
or later than 

Indeed, if exercise is delayed, then t (which is the “present” in the new 
pricing scenario) may be greater than i; in Q . In that case Q makes no sense. 

However, with appropriate modification the pricing argument may be adapted 
for t > 0. At time t > 0 the estimated exercise date is the value of s (f < s < r) 
which maximises 


gt(s) := E[pt (s, Zt)] = [ ptis,ZT)F{I[N]). 

If the distribution function F for the share price process Z is , 


= L 


p-(s-*) 


{z{s)-K)g^'^{i[N])- 


( 8 ) 

becomes 

( 9 ) 


As before, since the point-integrand pt is cylindrical, this gives qt{s) 

( f — du- 

\J-oo cry 27r(s — t) 




— K 

With ^ = 0, CT = 0.5, K = l.I, r = 2, and t = 1.5, Figureis the graph of qt{s) 
from s = 1.5 = t to s = 2 = t. This option has the same parameters (/r, a, k, t) 
as the one described in Figures and above, except that it is entered into at 
time t = 1.5 instead of time 0. But whereas Figure indicates an estimated 
exercise date « 1.83, Figure [^indicates an exercise date <Ji ,5 = 2. 


or 


dtis) = 


— p-(s-b / pM(s-t)+ 
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Thus the estimated value w(t) (or wt) of an American call option at time t 
{0 < t < t ) is Wt = 

In ^ + (p + 1^2) ^ / In t + (p - - t) 

cry/<;t - t j y (Ty/(,t - t 

(10) 

At time t = t the value of an unexercised American call option is 

Wr = maxjzT- — K, 0}. 

The preceding investigation uses distribution function . But the assumption 
that price processes Z follow a geometric Brownian distribution is generally 
false, as demonstrated in [2] Section 9.9, pages 479-485. 

Also, both European and American options have exercise values which de¬ 
pend on the value taken by the underlying asset on the date of exercise of the 
option. In contrast, options of the Asian type have value which depends, not 
just on the share value on the option exercise date, but on the average value of 
the share over the lifetime of the option: 

w { t ) = max{A(z7’) — k, 0}, 

where A{zt) is the average of the known or observed values taken by the share 
in advance of the exercise date r, which, for this purpose, is the “present”. 

So if the lifetime of the option consists of one million “time-ticks”, by time 
r the share value z{s) will have been observed for each of one million values of 
z(s), and, with T =]0,r], A{zt) is the sum of a million values of ^(s) divided 
by a million. 

However, at any time t earlier than r, including time t = 0 when the op¬ 
tion contract is created, the future values z(s) (0 < s < r) are unknown and 
unpredictable. But if their joint distribution F is known, then A{Zt) is a joint- 
contingent observable, 

A{Zt)c^A{zt) [R^,F] , 

and likewise friZr) — inax{A{ZT) — k, 0} F]. The values Wt of the Asian 
option at times t, 0 < t < t can then be deduced from (8.47) of [2], page 440: 

Wt = [fr {Zt)] = [ /,(zt)F(/[1V]). 

The defects of the hypotheses (assumption of geometric Brownian motion) un¬ 
derlying the preceding argument are outlined in Section 9.9 (pages 479-485) of 
[M] . Sections 9.10 and 9.11 of [2] outline an empirical method for estimating 
the actual joint distribution functions for such processes. 

The remainder of this article applies this method to estimate the value of 
an Asian call option on Glanbia shares, using the Maple program. 

Both European and American options have exercise values which depend on 
the value taken by the underlying asset on the date of exercise of the option. 







In contrast, other kinds of option contract have value which depends, not just 
on the share value on the option exercise date, but on the values taken by the 
underlying asset at various times during the lifetime of the option contract. 

For instance, an option of the Asian type depends on the average value of 
the share over the lifetime of the option: 

w{t) = max{^(z7’) — k, 0}, 

where A{zt) is the average of the known or observed values taken by the share 
in advance of the exercise date r. 

So if the lifetime of the option consists of one million “time-ticks”, when 
time T is reached the share value z{s) will have been observed for each of one 
million values of z{s), and, with T =]0, r], A{zt) is the sum of a million values 
of z{s) divided by a million. 

However, at any time t earlier than r, including time t = 0 when the op¬ 
tion contract is created, the future values z(s) (0 < s < r) are unknown and 
unpredictable. But if their joint distribution F is known, then A(Zt) is a joint- 
contingent observable, 

A{ZT)zzAizT) , 

and likewise fr{ZT) — max{A(ZT) — k, 0} [R!^, F] . The values Wt of the Asian 
option at times t, 0 < t < t can then be deduced from (8.47) of [T3], page 440: 


Wt = [fr (Zt)] = [ Mzt)F{I[N]), (11) 

in accordance with the risk-neutral martingale pricing argument of |14j Sec¬ 
tions 8.13-8.19, pages 433-466, with 


frizr) = max{A(zT) - k, 0} 


( 12 ) 


for Zt & R+. In the case of European and American options, fr depends on a 
single variable Xr or with r (or ij) fixed, so the latter integral reduces to an 
integral on R” (or R+). 


But the integral in (11) latter integral is infinite-dimensional. Note that if 
t > 0 is the “present”, some values of zt, = {- 2 ( 5 ) : 0 < s < r}, are known values, 
while the others are still in the future, and must be regarded as unpredictable 
observables. So 


rat = [ U {{zt>, zt>>)) F{I[N]), 

where the values zt' G are known and the values zt" £ are potential 

occurrences of the joint-basic observable 


Zt" — zt" 




^'^Unlike the assumption of log-normality, these arguments seem to be fairly reliable. 
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To sum up, 

wo = e~P^'"'> [ fr{zT)F{I[N]) where frizr) = max{A{zT) - k}. (13) 

Instead of assuming F = , pages 486-489 of [Hj suggest an empirical ap¬ 

proach using the risk-neutral, no-arbitrage pricing argument of |14j Section 8.14, 
pages 436-438. Applying this argument to the Glanbia share price series [224], 
the task can be tackled in steps such as the following. 

• Using the historic Glanbia data zt, construct empirical distribution values 
J^(/[iV]) as in Section 9.11 of [H] pages 487-489. 

• Estimate the growth rate /i of the price process. This can be done by, for 
example, calculating least squares regression for the Glanbia data zt- 

• Galculate a discounted version z'{s) = e~^^z{s) of the Glanbia price data. 

• Using the growth rate estimate /i, adjust the distribution values F so that, 

with respect to amended values F, the discounted prices z'{s),w'{s) are 
occurrences of martingales Z'j, ~ F], ~ [R[(], F]: 

E^[Z(] = z'(0), E^[IUi]=n;'(0) 

for all s, 0 < s < r; with z'{0) = z(0), w'{0) = w{0). 

• Galculate an estimated value rc(0) by estimating the integral 

n;(0)=n;'(0)=E^[IU;]= / n;'(r)F(/[7V]) = [ /.(4)F(/[1V]). 

For Asian options the latter integral depends on all occurrences z(s), 0 < s < r. 
So the integral is not cylindrical, and in order to estimate it, suitable regular 
partitionaH of = R[^, can be used, along with step function estimates of 

the integrand. 

Step function approximations are, in principle at least, easy enough to un¬ 
derstand and to apply—even in the case of multi-dimensional integrals, when 
we do not have recourse to the immense simplification that occurs when the 
integrand is a cylinder function whose integral happens to reduce to a one¬ 
dimensional integral, as in the cases of European and American options. 

But, even when step function simplicity is present, it is still necessary to 
engage with the financial principles involved in risk-neutral valuation. This 
confronts us with plenty of tricky issues. 

The main issue is converting the process distribution F to a martingale 
distribution F. How can this be tackled in the case of the Glanbia price process, 
for instance? 

(r, g)-partitions are a bit cumbersome in practice, and other kinds of regular partitioning 
will be used below. 
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Accepting that F cannot be taken to be the geometric Brownian function 
m Sections 9.10 and 9.11 (pages 486-489) suggests using a counting or 
“relative frequency” method to estimate F for an observable Xt — F], 

This is problematic. Suppose an observable X has distribution function P 
on sample space R, X ~ a;[R, P], so a P-measurable set of potential occurrences 
A C R has probability P{A). The relative frequency idea implies that “multi¬ 
ple occurrences” of A, expressed as a fraction of total occurrences, provide an 
estimate of P{A). Some justification of this estimate is provided by the Law of 
Large Numbers ([14] Section 5.9, pages 224-233). 

The trouble with this is that “multiple occurrences” of a random variable X 
(or of a process Zt) are contrary to the usual mathematical understanding of 
a random variable, or of a process. A random variable/observable is generally 
understood as a mathematical description of a single future possible occurrence 
of something; that is, a single future observation or measurement which has 
multiple possible values, just one of which can actually “occur”. 

So there are no “multiple occurrences” from which relative frequency can be 
calculated as estimate of probability. 

On the other hand, multiple occurrences, and relative frequency, can be 
somehow extracted from the Law of Large Numbers by postulating, not a single 
random variable with multiple occurrences, but multiple random variables, each 
with a single occurrence. 

But in order to connect this line of reasoning to the “relative frequency 
approximates to probability” idea, this family of random variables must consist 
of statistically independent random variables. 

Glanbia (or any other) historic share price data constitute a single occurrence 
Zt of a family of joint observables {Zt : t G T}, written as process Zt — 
ZT [R^,F]. The distribution function F(/[A^]) for the process can be thought of 
as consisting of transition probabilities for joint occurrences 


G I(tj), I < j <n), N = [t 


l5 ■ • • 5 


tu}- 


Can relative frequency be used to estimate F{I[N])1 This question is addressed 
in m pages 486-489. The idea is as follows. 

The Glanbia share price series begins with the end-of-day price 0.8 on Friday 
8 March 1991, and continues with the daily price for each day of 5-day weeks 
until Wednesday 9 March 2011 when the end-of-day price was 3.692. Thus the 
series contains over 5,000 individual items. 


The price transitions in /[A^] of (11) are possible transitions in the future. 
Assume that, at any stage in the “history” (past, present, or future) of the data, 
the likelihood of such a pattern of transitions is some more-or-less hxed amount 
for this particular series of share prices. 

It is possible to count how many transitions in the historic or past data 
correspond to the set of possible future transitions represented by I[N]. To 
convert this to a relative frequency, divide this number by the total number 
of possible transitions corresponding to I[N] in the past data. Provided the 
number of price transitions n in the event I[N] is small compared to 5,000, the 


11 



preceding assumptions suggesl[^ that this relative frequency could be used as 
an approximation to 

Next, to perform risk-neutral valuation, the (estimated) distribution function 
values J^(/[-/V]) must be converted to (estimated) martingale distribution values 
-F'(/[N']) for I[N] G I(R^). How can this be done? 

Consider a single share price Zt at some fixed time t, 0 < t < t. For cells 
It =]M,t’] G I(R+), methods such as those described above can be used to find 
an estimated distribution function Ft{It) for the observable Zt ~ Zt [R+, Fi], 

V {zt G It) = Ft{It), ^jtt=F[Zt]= [ ZtFt{It) = 

The basic observable Zt can be regarded as a contingent observable in R^, with 
/(2't) = Zt, /{Zt) ^ f{zT)[Rl,F]. Then / is a cylinder function, /{zt) = Zt, 
t fixed; and 

E[Zt] = I f{zT)F{I[N]) = [ ztFtilt) = fit 

Since F{I[N]) = Ft{It) for N = {t}, /[iV] = If 

Given a price process Zt — zt [R+ , F], risk-neutral pricing theory requires 
the construction of a process Zt — z^lR^, F] such that, under the distribution 
function F, the discounted values e“'’®z(s) are occurrences of a martingale, with 

[e-P^Zs] = z{0), or [Zs] = e'’"z(0) 

for each s G T, where z(0) is the known initial price, and where, for present 
purposes, the risk-free interest rate p is taken to be constant for all s. 

Example 1 The construction of Ft can be illustrated as follows. Suppose we 
have an observable X ~ x[R, i?], with mean value p, = E[X]. Suppose we wish, 
by changing the values of H, to construct a different observable X whose mean 
is some given number E^[X] = p. This can be done in two stages. First 
construct X' ~ x'[R, FT] with mean 0, then construct X = a;[R, fl] mean 

p. For the first stage, suppose I' = ]u',v'] G R. Write 

u = u' + p, V = v' + p, so u' = u — p, v' = V — p, 

and, for It =]u,v], define Hfl') = II{I), so E^ [X'] = 0. The second stage is 
similar. For It = ]u,v] G R write 



u = u + p — p, V = V + p — p, so u = u — p + p, V = V — p + p; 
and, for It =]m, u], define H{I) = F(/), so E^[X] = p, as required. 

^^This can be criticized from various angles. One benefit, however, is that there is no 
presumption of independence. Any conditioning or dependence is picked up by the counting 
process. 

^^The elements x,x',x are not different numbers; each of them represents an arbitrary 
possible occurrence in R,. 
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Returning to the construction of Zt — zr[R+,-F], Example ^ shows how to 
construct Ft for each t GT, so, for N = {t} the values F(J[A^]jcan be found. 

But if Zt is to be defined as a process or joint observable, the joint distri¬ 
bution function .F'(/[iV]) must be defined, not just for singletons N = {t}, but 
for all N G Af{T) and all corresponding I[N] G I(R!^). The counting procedure 
addresses by considering various patterns of transition. 

The Glanbia daily price share data described in El Section 9.11 (pages 487- 
489) consists of 5225 daily end-of-day Glanbia share prices, beginning Friday 8 
March 1991 and concluding Wednesday 9 March 2011. 

Suppose the latter point in time is the present, and a 3 month call option 
contract on Glanbia shares is entered into at that moment. Since the share 
prices are reported daily on a 5-day week basis, take the 3-month term of the 
option to be 60 days, starting on morning of 10 March 2011 and expiring at 
end-of-day Wednesday 1 June 2011, a total of 5219 individual days/prices. 

Suppose the option is the Asian type, with pay-off depending on the average 
share price over the term of the option. 

To find an economic price for entering into this contract, it will be useful to 
have estimates of the daily expected value of a share during the 60-day term of 
the option; likewise the daily standard deviation of the share price. Also, joint 
probability values E(/[iV]) need to be estimated for the share price process 
during the term of the option. 

Maple calculations on the historic share prices can be used to make such 
estimates. 

The historic Glanbia price data [4] can be found in an Excel file in 
https://sites.google.com/site/stieltjescomplete/ 

The immediate aim is to produce estimates of expected share price value 
for each day of the term of the option—that is, the expected share price on 10 
March 2011, on 11 March 2011, and so on, up to and including 1 June 2011; all 
of which are in the “future”. 

Likewise, to estimate daily share price standard deviation for each day of 
the term of the option. Maple calculations on the historic or past data provides 
such estimates. 

Calculation 1 

1. with(ExcelTools); with(Statistics); 

2. Q Import("glanbia.xls") : 

3. A1 ■- seq(g[j][2],J = 2500..5010) : A2 := [Al] : A3 :=Vector(A2) : 

4. J1 :=seq(j,j = 2500..5010); J2 ;= [71]; 73 ;=Vector(72); 

5. y :=ExponentialFit(73, A3, t); 

6. GlanbiaGraph ;=[seq([t, Q[t + 7][2]], t = 2500..5000)]; 

7. plot ( [GlanbiaGraph, y],t~ 2500..5000) 
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8. p := 7;q ■— p + 59; 

9. for k from 1 to 5160 do 

10. A[k] ~ seq(g[j][2],j = p..q)\ 

11. B[k] ■.= [A[k]\\p--p+l\q--q+l\ 

12. m[k] Mean(i3[fc]); s[fc] StandardDeviation(i3[fc]): 

13. end do; 

14. m := seq(m[fc],fc = 1..5160);s :=seq(s[fc],fe = 1..5160); 

15. S ;= seq(s[A;],fc = 5000..5160);T ;= [S'];s0 ;= Mean(r) 

16. M[l] := 3.7*exp(0.006);M[2] := 3.7*exp(0.012); M[3] := 3.7*exp(0.018): 

17. for k from 1 to 5219 do 

18. z[k] := Q[fe + 6][2] : 

19. end do: 

20. N := 6; J := seq{j,j = —3..2); seq(g, i = 1..3); seq(b, i = 1..3); 

21. for i from 1 to 216 do 

22. for j from 1 to 3 do 

23. q[j] := trunc((i - l)/iV^i - 1)) + 1; 

24. b[j] ■- mod((7[j] - 1, A^) + 1; 

25. p[i\[j] := J[b[j]] 

26. end do 

27. end do 

28. ex := 3.5; wO 0; 

29. Prob ;= 0; 

30. for i from 1 to 216 do 

31. /:=0; 

32. for k from 4550 to 5000 do 

33. if {z[k + 20] >= m[k + 20] + s[fc + 20] * p[i] ]1] 

34. and z{k + 20] < m[k + 20] + (p[*][l] + 1) * s[fc + 20]) 

35. and z\k + 40] >= m[k + 40] + s[A: + 40] *p]i][2]) 

36. and z[k + 40] < m[k + 40] + (p[i][2] + 1) * s[k + 40]) 
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37. and z[k + 60] >= m[k + 60] +p[*][3] * s[l + 60] 

38. and z[l + 60] < m[l + 60] + (p[i][3] + 1) * s[l + 60]) 

39. then /;=/+!: 

40. end if; 

41. end do; 

42. P:=(l/50)*/; 

43. Prob ;= Prob+P; 

44. al ;= (M[l] +p[i][l] * sO + M[2] +p[i][2] * sO + M[3] +p[i][3] * s0)/3 — ex', 

45. if al > 0 then a2 ;= al; end if; 

46. if al <= 0 then a2 := 0; end if; 

47. a3 := a2*exp(—0.02); 

48. w := a3* P; 

49. wO ~ wO + w: 

50. end do; 

51. wO; 

52. evalf(Prob); 

For the purpose of explanation, the lines of Maple code above are numbered. 
These numbers are not part of the program and should not be included in actual 
Maple code. 

Lines 1-7 take the historic Glanbia data and import them into the program 
for the purpose of performing Maple calculations on them. These data can be 
found in an Excel file in https://sites.google.com/site/stieltjescomplete/ 
and should be stored in the same folder in which the above Maple code is located, 
in order for this code to be able to access them. In version 18 of Maple, the file 
with Import command must, on first use, be closed and then re-opened before 
Import will work successfully. 

Lines 3 and 4 convert the share price data to a format suitable for line 5, 
using about half of the available prices—from day 2500 to day 5010; it is a 
matter of judgement which data to include. 

Line 5 calculates a least squares regression for the price data, in the form 
y = a exp (/it). The idea here is that the prices follow some underling growth rate 
/i of a “proportional” character; with random variation of prices superimposed 
on this underlying trend. In other words, if the share price on day t is Zt, and 
disregarding the superimposed random variation of prices, then 

zt+i . . ^ f zt+i \ 

- IS constant, in - = u. 

Zt \ Zt J 
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In this case the exponential best fit calculated by Maple is approximately 

y = 

so daily growth rate is approximately /r = 0.0007. This translates into an annual 
growth rate of the share values amounting to about 20%. 

Lines 7 and 8 of the Maple code produce a graph {y) of this underling growth, 
superimposed on a graph of the actual prices. See Figure at the end of this 
article. 

For the purpose of pricing an Asian option entered into ’’today” (9 March 
2011), when “today’s” share price is 3.7, a risk-free daily growth rate of 0.0003 
(which is less than the underlying growth rate /r = 0.0007), so the underlying 
trencp^ values of the share during the period of the Asian option (up to 1 June 
2011, or 60 days in total) is taken to be 

y = 0.1e°°°°^*. 


Then, provided we can superimpose on these price trend values the appropriate 
amount of random variability, we can estimate the value of the Asian option by 
using (11). 

Lines 8-13 are the first step in accomplishing this. This part of the program 
calculates successive 60-day average values of the historic share prices. Also 
the corresponding standard deviations. The latter demonstrate the scale of the 
random variability or volatility in these share prices, and are a kind of indicator 
of the joint likelihood distribution F{I[N] that “determines” the share price 
process Zx — 2 :t[R^,F’]. 

These “moving averages” and “moving standard deviations” are illustrated 
in Figures and at the end of this article. 

The Maple program will also use these mean-and-standard-deviation data 
as “3 (t” partition points in a regular partition of , T being the 60-day period 
or term of the Asian option, in order to estimate the option valuation (11). 

Line 15 takes an average sO of the 60-day standard deviations. This is to 
be taken as the standard deviation of the daily share price for each day of the 
60-day term of the Asian option—see line 44. 

Line 16 computes three 20-day trend values M[l],M[2],M[3] for the 60- 
day term of the option. But these are not “true” trend values; they are the 
hypothetical values showing a risk-free daily growth rate of 0.0003 needed to 
carry out the risk-neutral martingale valuation of 

Lines 17 to 19 ascribe the familiar notation z[k] (or Zt) to the historic share 
price values. Lines 20-27 compute the permutations, with repetition, of the 
six numbers —3, —2, —1, 0,1,2 taken three at a time. There are 6^ = 216 such 
permutations, which, for i = 1 to 216, are each given by 


{pim, pimi pwp]). 

^^Not the actual trend, just the hypothetical risk-free trend needed for martingale pricing 
of discounted values. Line 16 of the code shows these successive “risk-free trend” values at 
intervals of 20 days during the term of the option. 
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One such permutation is (1, —1, 2); that is 

p\i][l] = l, p[i\[2] = -l, p[i][3] = 2. 

With mean m and standard deviation s, these numbers are then used to con¬ 
struct 3cr partitioning intervals 

[m - p[i][j] *s,m- {p\i][j] -b 1) * s[, 

which, for p[i\[j] = —3, —2,..., 2, give six cells 


[to— 3s, m—2s[, [to— 2s, TO—s[, [to— s, to[, [to, to-|-s[, [to-|-s, to-|-2s[, [to-|-2s, to-|-3s[. 

On any particular day the share price will usually be found to lie in one of 
these intervals. The probability of \zt — to| > 3s is small, and can be ignored 
in the option valuation calculation. Take the start of 9 March 2011 to be time 
t — 0, and 1 June 2011 to be t = t, and T =]0, t], with t = n = 20, t = T 2 = 40, 
and t = r = 60. For j = 1, 2,3, write 


= Wij] + p[i] [j] * sO, M[j] + {p[i] [j] -b 1) * s0[ 


(14) 


where the values sO and M[j] are given by lines 15 and 16 of the Maple code. 

Then, with N = {ti,T2,t} the domain R[([ of the option is partitioned by 
cylinders 

( 15 ) 

1=1 

There are 216 such cylindrical cells, corresponding to the 216 permutations 
{p[i] [1], p[i] [2], p[i] [3]). This is a regular partition in R[^. It does not fully exhaust 
the domain R[(), but the probability that zt is not in one of these cylindrical 
intervals is very small. Line 28 of the Maple program sets the exercise price of 
the option at ex = 3.5. The assignment wO := 0 sets the initial value of the 
Riemann sum estimate of (11). The Riemann sum has 216 terms, one for each 
permutation (p[i][j];j = 1,2,3). Each term of the Riemann sum has value w 
(line 48), and the term values w are accumulated as wO in the 216 iterations of 
lines 48 and 49. 

The ultimate value returned by wO is the value of the Riemann sum estimate 


of (111. For option exercise price ex = 3.5, the Maple program gives initial value 


wO = 1.31. This is the Wg or r(;(0) of (111. The expression Prob of line 29 is 
used in the code to keep track of, or accumulate, the probabilities generated 
in the calculation. Total probability returned by line 43 is approximately 0.97, 
indicating that the 3cr regular partition of R[^, though non-exhaustive, is almost 
full. 

Using the historic Glanbia share price data, lines 30-43 of the program use a 
“relative frequency”-type argument to produce estimates of the joint transition 
probabilities F{I[N] needed to calculate the Riemann sum estimates wO of 0- 
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Line 32 shows that this program only uses 50 cycles or iterations of the 


3 X 20-day transitions in the partitioning cells (15). For greater accuracy, it is 
easy to increase this sample size. 

The relative frequency calculation for a single permutation of transitions 
(p[*]b] ■ j = 2, 3) is done in line 42. This value of P is the estimate of 

for a single term of the Riemann sum for Q- 

For i = 1,..., 216, each permutation (p[i][j] : j = 1, 2,3) generates a single 
term of the Riemann sum for The calculation of each such term is done 

in lines 44-48. 

How are the values F"(/[iV]) converted to risk-neutral probabilities F{I[N])? 
This is accomplished in line 44. The quantities M[j] -|-p[z][j] * sO are the lower 
bounds of partitioning cells of the form (14) and ( |I^ , and these are taken as 
tag points or evaluation points for the function A{zt) in 0- 

But these boundary points are unlike the corresponding boundaries in lines 
33-38. The latter are growing at the daily trend growth rate = 0.0007, whereas 
the bounds in line 44 are growing at the risk-free growth rate p = 0.0003. 

So, in some sense, the probabilities F"(/[iV]) produced from lines 33-42 are, 
in line 44 applied to the “wrong” cells of . Except that this is the appropriate 
adjustment needed to produce the hypothetical F(J[A^]) in R!^. 

The point-integrand in ( [TT|) is the function ( |l^ , and this is calculated as a2 in 
lines 44-46. The discounted^ value of (12) is calculated in line 47. Multiplying 


this by P = P(/[iV]) gives the value of a single term w of the Riemann sum. 
Summing over all the permutations, i = 1,..., 216, gives the Riemann sum 
estimate of the option value at initial time t = 0: wO = wq = w(0) = 1.3. 

To test out the above Maple code, it is advisable to test its component parts 
individually in order to more easily detect and correct transcription and coding 
errors. Also, at the relevant points in the code, alternative parameters and 
estimation tactics can be experimented with. 

In the above code the average A{zt) was estimated at only three 20-day 
intervals. For a 60-day option the Glanbia share price data in the Excel file 
permit up to 60 daily price values to be used to calculate A{zt)- In practice, 
a reasonable balance must be struck between the demands of accuracy and the 
scale of computing power available. The latter quickly escalates when more 
accuracy is demanded. 

[14], page 480, has a twenty-year graph of the Glanbia share prices. Figure]^ 
is a graph of part of the data, with superimposed exponential regression graph 
y = exp(0.0007t). 

The first average calculated by Lines 8-13 of the Maple program is the mean 
of the prices for day 1 (8 March 1991) up to and including day 60. The second 
average is for day 2 (9 March 1991) up to and including day 61 of the data. 
The final average is the mean of the prices on day 5160 of the data up to and 
including day 5219 (9 March 2011). 


^®For 60-day discounting, the daily rate 0.0003 becomes approximately 60 X 0.0003 = 0.018. 
Line 47 uses instead of for discounting at the risk-free rate. 
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Figure 7: All moving averages Figure 9: All moving SB’s 




Figure 8: Exponential regression Figure 10: Option values 


This is a total of 5160 cycles, giving 5160 means and 5160 standard devia¬ 
tions. The latter are displayed in Figure]^ 

Since the first moving average is the average for the first 60 days, it can be 
applied to day 30 or day 31, half-way through the first cycle. Likewise for each 
of the other moving averages; and also the standard deviations. 

The graph of moving averages. Figure is a “smoothed out” version of 
the graph (Figure]^ of the original Glanbia share price data. But it is not as 
smooth as the exponential regression graph for the data, in Figure 

Figure [T0| shows how the initial value of the Asian option depends on the 
terminal exercise price. A low exercise price of 3.0 results in high option price 
of about 5.5. A high exercise price of 3.8 produces a low option price of about 
0 . 02 . 
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